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Topological charge distributions in 2 dimensional CP^ model with S-term is calculated. 
In strong coupling regions, topological charge distribution is approximately given by Gaus- 
sian form as a function of topological charge and this behavior leads to the first order phase 
transition at = tt. In weak coupling regions it shows non-Gaussian distribution and the 
first order phase transition disappears. Free energy as a function of 6 shows "flattening" 
behavior at 9 = Of < tt, when we calculate the free energy directly from topological charge 
distribution. Possible origin of this flattening phenomena is prensented. 



§1. Introduction 

Quantum Chromodynamics(QCD) admits the topological term, i.e., ^-term. It 
leads to "strong CP violation" . Experimentally, the magnitude of 9 parameter which 
controls the magnitude of strong CP violation is severely limited to a tiny value ( |^| < 
10~^). Lattice gauge theory( LGT), formulated by K.G. WilsonB) and developed by 
M.CreutzB', has made a realistic progress in understanding the nature of QCD. 
In LGT in which Euclidean space-time is adopted, usual gauge coupling is purely 
real quantity. On the other hand, O-teim appears as a purely imaginary quantity. 
Combining both of these leads to complex numbered coupling. The action given 
by this complex coupling constant gives complex valued Boltzmann weight. Except 
for a few number of cases, the works in LGT are done for real couplings. It is, 
however, very important to understand the role of complex coupling cases, i.e., the 
system with O-teim. Monte Carlo simulations are confronted with a difficulty due to 
complex Boltzmann weight. Monte Carlo simulation in this case was made possible 
to some extent by Bhanot et al.&i), Wiese& and Karliner et al.i'B. 

The method is summarized as follows; l)First we obtain the topological chargei^ 
distribution P{Q) with the use of the Boltzmann weight defined by the real coupling 
constant. Fourier series 

Z{e) = Y,P{Q)e'^Q 

Q 

gives the partition function Z{9) for the system with 0-term. 
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2) There are some technical problems to obtain topological charge distribution. 
Since the topological charge distribution is a rapidly decreasing function, we use (i) 
"set method", in which the whole range of Q is divided into number of sets. In 
each set, MC simulation is performed, (ii)even after dividing into sets, the distribu- 
tion is still rapidly decreasing, so we use "trial distribution method", in which the 
Boltzmann weight is normalized by some appropriate trial function so as to realize 
the slowly changing topological charge distribution. This method was applied to 
2 dimensional U(l) system with 0-term . Wiese showed that the first order phase 
transition is observed at = vr for all real couplings B^. 

Two dimensional CP^^^ system has many qualitative similarity with 4 dimen- 
sional QCD, like, asymptotic freedom, confinement, deconfinement phase( Higgs 
phase) and topological excitations( instanton). 

Schierholz studied 2 dimensional CP^ system with 0-terma). He observed that, 

1) in strong couplings, the first order phase transition occurs at a critical value 
6c = TT and that 

2) in weak couplings, deviation of the position of the first order transition from 
6c = TT to smaller values. He conjectured that 6'c ^ as /3 ^ oo( weak coupling 
limit) and that the continuum limit value of 6 is zero. 

If this conjecture is valid, it gives quite important physical result, but almost 
no confirmation is made by other groups. One of the aim of the present paper is to 
investigate this issue 0^0. 

We will briefly summarize the approach other than MC. We analyzed 2 di- 
mensional U(l) system with 0-termtlP based on the group character expansion 
method Topological charge distribution P{Q) is shown to be given by a Gaussian 
function exp(— KyQ'^)) foi' ^ It leads to the result that the partition function is 
given by the third elliptic theta function 'd^{v, r). It has a remarkable property that 
the partition function has an infinite numbers of zeros as a function of z = e*^, where 
6 is extended to complex values. In infinite volume limit, these infinite numbers of 
partition function zeros accumulate to z = — 1 according to Ky oc 1/V. Namely, 
first order phase transition is expected at = vr in the limit V ^ co, where V = LP' 
denotes the volume of system. 

Foiir dimensional Zi\j gauge system with 0-term is discussed by Cardy and Rabi- 
noviciB. They pointed out that ^-term leads to symmetries under dual transforma- 
tion of Hamiltonian and this symmetry leads to very rich phase structures, namely, 
confinement phase, Higgs phase( deconfinem ent phase) and condensation of electric 
and magnetic charges( oblique confinement) US. Their argument is based on free 
energy and rather qualitative. It will be of much interest to investigate this system 
by MC simulation or other method( e.g., renormalization group method). 

In the previous paper 113, we performed MC simulation on CP^ system with 
0-term in strong and weak coupling regions. In strong coupling regions, first order 
phase transition is observed at = 7rll§'©. In weak coupling regionstl^ we observed 



*' After completion of the main part of this paper, we are informed that similar investigation 
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that topological charge distribution deviates apparently from Gaussian form. Free 
energy in weak coupling regions is approximately given by cos^ form, which shows 
only a few small Q's control the topological charge distribution. In intermediate 
coupling regions, obtained results were not so clear. In order to avoid the effect of 
statistical error as far as possible, we used fitted analytic form for P{Q) to obtain 
Z(9),F(9) etc. Then even in weak coupling region we observed no flattening behavior 
but change of behavior of Z(6), F{9). The power series fit to P{Q) has the property; 

1) for /3 small( strong coupling regions), it is simply given by Gaussian form( 
exp(-Ky(52)). 

2) for /5 large( weak coupling regions), we need term in the polynomial fit of 
the exponent of P{Q). 

3) for (5 intermediate, power series fit suffers from quite large value. 

In this paper, CP^ model with 0-term will be studied by MC methodic. The 
emphasis is put on intermediate coupling regions(/3 = 3 ~ 4). From P{Q) we obtain 
Z[9) and we obtain free energy per unit volume by F{6) = —{\/V)\n.Z{6). F{9) 
grows smoothly in smaller 9 and shows "fiattening" at some value 0/, beyond which 
F{9) becomes almost fiat. For example, at /? = 3.5 and L = 20( 100000=100 kiro 
iterations), fiattening is observed at 9 = 9f ^ O.STtt. Does it correspond to the "first 
order phase transition" found by Schierholz in GP^ system? 

Does the fiattening mentioned above really mean a phase transition? We should 
be careful about the statistical error induced by the simulation. For example, after 
N times of sweeps( =iterations), we would have statistical error of order 

5P = |error| ~ l/\/iV. 

Note that 5P is defined as positive. 

Set method and trial function method will give us P{Q) containing great number 
of orders of decrease from Q = to large value of Q. If observed P{Q) is fitted by a 
smooth function Ps{Q), it can be expressed as, 

P{Q) = Ps{Q) + ^P{Q), 

where AP{Q) is the difference between observed P{Q) and the fitted function Ps{Q)- 
Quantity |Z\P((5)| will be approximately given by 5P x P{Q). The quantity AP{Q) 
can be both positive and negative( or zero), since it is a fluctuation. The value of 
|Z\P((5)| is approximately of order of {\/^fN)P[Q) and AP{Q) itself is expected to 
be small at large Q. Partition function Z{9), however, is given by the fourier series 
and the error in evaluation of Z[9) is controlled by the largest one, i.e., |zAP(0)| ~ 
5P ~ 1/yfN . Topological charge distribution P{Q) is a rapidly decreasing function. 
When we estimate Z{9) from 

Z{9)=Y,P{Q)e'"^, 
Q 

the value of Z{9) will also be a rapidly decreasing function of 9 and the quantity 
Z\P(0) will play an important role at large 9. Since the partition function obtained 
by Ps{Q) is a rapidly decreasing function, true Zs{9) will be masked by the error 
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AP{0) at large 9, which is much larger than PsiQ) itself at large Q . This effect 
leads to the shape of Z{9) with flattening at some Of just like the illustration shown 
in Fig.l{6f ~ O.SSvr in this figure). 



' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 

0.5 1 1.5 2 2.5 3 3.5 
theta 

Fig. 1. An illustration of partition function as a function of theta. Direct fourier series( curve with 
flattening behavior) and one( Gaussian curve) calculated from smooth fit to P{Q). 

The value of Z{6) will almost coincide with Zs{0) m. < Of but it deviates far 
from Z{0)s va. > Of. Zs{0) will become quite small as ^ — > vr but obtained Z{6) 
does not show such suppression. 

In this way the combined method of set and trial function leads to satisfactory 
result for P{Q) but we still encounter difficulty in obtaining correct Z{0) at large 
0's. 

In order to confirm the physical flattening (if it exists) , we need enough accuracy 
of Z{0) over huge range of magnitude. For example to confirm the behavior of 
partition function down to 10^^, we need ~ 10^^ iterations. 

To obtain information about flattening within available computer facilities, we 
should analyze at larger (3 or smaller volume V , since P{Q) decreases faster at /? = 
large or at y = small, partition function Z{0)^ which is a Fourier transform of P{Q), 
falls slowly at /3 = large orV= small. In these regions Z(9) will be reliably estimated 
over whole region of ^(0 < ^ < vr). 

§2. CP^ model in two dimensions 



2.1. CP^ model in two dimensions 

We define the CP^ model with 0-term. in the two dimensional euclidean lattice 
space. The CP^ model consists of three component complex scalar fields, 

/ zi{n) 



Z2{n) 



(2-1) 
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at each site n, and constrained by 

3 



5^<(nK(n) = l. (2-2) 



a=l 

Action, which is U(l) gauge invariant, is given by 

Se[z] = S[z\-ieQ[z\ (2-3) 

^ = /5EE{l-E l<H^a(n + (2-4) 

n=l /i=l,2 a=l 

where Q is topological charge, and ^a(n) is coupled with Zaiji + /u) at the nearest 
neighbor site n + /x(jLt = 1,2). Topological charge Q is defined by 

Q = V e^^(yl^(n) + A^(n + ii) - A^{n + u)- Ay{n)) (2-5) 

where, ^,i(n) = arg(2;*(ra)2;(n + /x)). Topological charge Q is given by the winding 
number, i.e., how often the field ^^(n) covers U{1) space when we move around the 
periodic boundary of the whole space time once. 

The partition function as a function of the coupling constant (3 and is 

^^^^ - !Dz*Dzc;^^{-S[z\) 
where Dz denotes integration over all fields, and the free energy F(9) is given by 

F{e) = ~\nz{e), (2-7) 

where V is volume of the lattice space. 

We would like to compute 9 vacua effects, but the complex Boltzmann factor 
e^^Q prevents us from performing ordinary Monte Carlo simulation. In order to get 

over this difficulty, we follow Wiese's idea. The updating is restricted to the fields 
in the topological sector. Thus the phase factor e^^^^^^ is replaced by a constant and 
can be factored out from the functional integral. We can perform ordinary Monte 
Carlo simulation using the real action S\z\. The partition function is given by the 
summation of the topological charge distribution P{Q) weighted by e*^'^ in each Q 
sector. In practice, we get the topological charge distribution P{Q) at ^ = by 
the Monte Carlo simulation. We obtain the partition function Z{6) by taking the 
Fourier transform of P{Q). 

Z{9) = J2e'^'^PiQ) (2-8) 
Q 

where P(Q) is 

^ jDz<Q)Dz(Q)eM-S[z]) 
jDz*Dzexp{-S[z]) 
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The integration measure Dz^'^^ is restricted to the fields in the topological sector 
labeled by the topological charge Q. Note that J2q P{Q) = 1- 

The expectation value of an observable O is given in terms of P{Q) as 

where {0)q is the expectation value of O at = for a given Q sector 

_ jDz^Q^Dz^Q^ O e-^M 
^ ~ J Dz<Q) DziQ) e-SM • ^ ^ 

2.2. Algorithm 

We measure the topological charge distribution P{Q) by Monte Carlo simulation 
with the Boltzmann weight exp(— 5), where S is defined by eq.(2.4). The standard 

Metropolis method is used to update configurations. To calculate P{Q), we count 
the number of times the configuration of Q is visited by histogram method. The 
distribution P{Q) damps very rapidly as | Q | becomes large. We need to calculate 
P{Q) at as larger \Q\ as possible, which would contribute to F{9),{Q)e and {Q'^)b 
because they are obtained by Fourier transformation of P{Q) and the derivatives of 
the partition function by ^. In order to obtain P{Q) at larger Q's, we apply following 
two techniques, 

a) the set method, 

b) the trial distribution method. 



The set method ; 

Range of whole Q is divided into number of sets 5, ( i = 1, 2, • • •). Monte Carlo 
updatings arc done in each set Si (Si = {Q \ 3{n — 1) < Q < 3n}). In the process, 
we start from a configuration within the set Si, and we produce a tentative configu- 
ration Cj. When Q of this tentative configuration Ct stays within one of the bins in 
Si, the configuration Ct is accepted, and the count of the corresponding Q value is 
increased by one. On the other hands, when Ct goes out of the set Si, Ct is rejected, 
and the count of the Q value of the old configuration is increased by one. This is 
done for all sets Si{ i = 1,2,- ■ •). 

The trial distribution method; 
Topological charge distribution P{Q) is still rather rapidly decreasing function of Q 
even in each set, the number of count in each set is sharply peaked at smallest Q and 
almost no counts in the larger Q's, so statistical weight is modified by introducing 
trial distributions Pt{Q) for each set. Namely Boltzmann weight is replaced by 
exp[— 5] /Pt{Q). This is to remedy P{Q) which falls too rapidly even within a sets 
in some cases. We make the counts at Q = 3(i — 1), 3(i — 1) + 1, 3(i — 1) + 2 and 3i 
in each set Si become almost the same. As the trial distributions Pt{Q), we apply 
the form 

Pt(Q) = Aexp[-^Q2] (2-12) 
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where the value of Ci{(3) and Ai depend on the set Si. That is, the action during 
updating is modified to an effective one such as ^eff = S + \n.Pt{Q). 

To reproduce a normahzed distribution P{Q) in the whole range of Q from the 
counts at each set, we make matching as follows: 

i) At each set Si{i = 1, 2 • • •), the number of counts is multiplied by Pt{Q) at each 
Q. We call the multiplied value Ni{Q), which is hopefully proportional to the 
desired topological charge distribution P{Q). 

ii) In order to match the values in two neighboring sets Si and we rescale 
Ni+i{Q) so that Ni+i{Q) Ni+i{Q) x r, where r = Ni{Q = 3i)/iVi+i(Q = 3i), 
the ratio of the number of counts at the right edge of Si to that at the left edge 
of Sj+i. These manipulations are performed over all the sets. 

iii) The rescaled SiS are normalized to obtain P{Q) such that 

P(Q) = . (2.13) 

i Q 

§3. Numerical result I 

We use square lattice with periodic boundary condition. We measure P{ff) 
in various lattice sizes {V = l?^ and coupling constants The error analysis 

is discussed in the previous paper©. To check the algorithm, we calculated the 
internal energy. It agrees with the analytical results of the strong and weak coupling 
expansions. Using the calculated -P(Q), we will estimate the free energy F{Q) and 
its derivative (Q)^, respectively. To obtain Z{Q\ we made two methods. 

(i) We use the measured Piff) directly to obtain Z{d^ ("direct method"). 

(ii) To avoid the error problems discussed in the introduction. We first fit the mea- 
sured P{ff) by the appropriate function -Ps(Q)) and then obtain Z{d^ by Fourier 
transforming from Ps((5)( "fitting method"). 

3.1. Topological charge distribution P{Q) 

In this subsection we discuss the topological charge distribution P{Q). In Fig. 
2, we show the measured P{Q) for various /3's for a fixed volume (L = 20). These 
P(Q)'s have a different behavior between the strong coupling regions and the weak 
coupling regions. In strong coupling regions, P{Q) exhibits Gaussian behavior. In 
weak coupling regions, P{Q) deviates gradually from the Gaussian form. In order 
to investigate the strong coupling regions in detail, we use the chi-square-fitting to 
In P{Q). Table I and II display the results of the fittings, i.e., the coefficients a„ of 
the polynomial J2n o-nQ^, for various /3's with the resulting x^/ dof, where dof means 
the degree of freedom. (i)For f3 < 1.5, the P(Q)'s are fitted well by the Gaussian 
form. (ii)For 2 < /3 < 3, terms up to the quartic one are needed for sufficiently good 
fitting. (iii)For 4 < /3, the fittings using the quartic polynomial turn out to be quite 
poor. Here we discuss the volume dependence. In strong coupling regions, P{Q) is 
fitted very well by Gaussian form for all values of V , 

P(Q) ocexp(-Ky(/3)Q2), (3.1) 
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where the coefhcient kv(/3)(= 02) depends on /? and V. Fig.^ shows lnKv{(3) vs. 
InV for a fixed /?(= 1.0). We see that Ky{(3) is clearly proportional to 1/V 

K = - (3-2) 

This 1/y-dependence of the Gaussian behavior determines the phase structure of 
the strong coupling regions. This result agrees with the CP^ result. Fig]^ shows the 




0. 001 I I I I I I I 

100 1000 roooo 

V 

Fig. 3. Ina2 vs. \nV for (3 = 1.0 and iV = 10^ 

volume dependence of P{Q) for L = 15, 18, 20, 25 at /3 = 3.45 . Parameters aj, 7 = 
2.0 at P = 3.45 (see Table. V) give F-dependence. We find slight deviation from 
exact l/V law, but a clear volume dependence is observed, which gives ~ 
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1-15 




1-18 




1-20 




1-25 



50 100 150 200 250 

Q"2 

Fig. 4. P(Q) vs. 0^ for P = 3.45 and iV = 5 * 10^ 

Table I. The results of chi-square-fitting to In P{Q) in terms of the polynomial ao + a2Q^ for various 
/3, where dof means degree of freedom. 





x' 


x7dof 


ao 


a2 





20.058 


0.872 


-2.509 


-0.021 


1 


83.936 


3.649 


-2.416 


-0.025 


1.5 


136.066 


5.916 


-2.368 


-0.028 


2 


769.635 


34.625 


-2.296 


-0.033 



Table II. The result of chi-square-fitting to In P{Q) in term of the polynomial ^ anQ" for various 
(3. 



p 


x' 


x7dof 


ao 


ai 


02 


as 


a4 


2 


27.835 


1.392 


-2.273 


0.006 


-0.035 


5.123 * 10" 


-b 


1.359 * 10" ' 


3 


30.561 


1.528 


-1.950 


0.01 


-0.069 


5.765 * 10" 


-4 


-2.904 * lO"*^ 


4 


3496.599 


174.830 


-0.791 


-0.684 


-0.188 


5.842 * 10" 


-3 


-8.643 * 10"^ 


5 


4027.247 


201.362 


4.147* 10"^ 


-4.468 


-0.159 


6.161 * 10" 


-3 


-1.015 * 10"* 


6 


5175.739 


258.787 


-3.341 * 10"^ 


-8.806 


-0.054 


1.423 * 10" 


-3 


-1.393 * 10"^ 



3.2. Free energy and expectation value of topological charge 

The partition function Z{6) is given as a function of 6 by (2-8) from P{Q)- The 
free energy is 

F(0) = -ilnZ(0). (3-3) 
The expectation value of topological charge is 
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In strong coupling regions, we have Gaussian behavior of P{Q), and the l/lZ-law 
appears to hold up to L = 70. It is natural to expect that this behavior persists 
to y — > oo. Let us look at how the 1/V-law affects F{9) and {Q)e- By putting 
C{(3) = 9.9, in P{Q) = exp[-^g2] p ^ ^ calculate F{e) and {Q)g from 

(2-8), (3-3) and (3-4). Figure.^ and ^ show their volume dependence. As V is 



0. 25 



.0. 15 




0. 05 



I " " I " " I " " I " " I " " I " " I " " I " " I 
0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 1 
theta/pi 





1 = 20 




1 = 10 




1 = 6 




1 = 4 



Fig. 5. Free energy vs theta in strong coupling regions. L = 4,6, 10, 20 



0. 16 




0. 1 0.2 0.3 0.4 0.5 0. 6 0.7 0. 8 0.9 1 



theta/p i 

Fig. 6. < Q > vs theta in strong coupling regions. L = 4,6, 10, 20 



increased, F{9) very rapidly approaches a quadratic from in 9 from below. Its first 
moment {Q)e develops a peak near 9 = n, and the position of the peak quickly 
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approaches n as V increases. The jump in {Q)e would arise at ^ = vr as y 
This indicates the existence of the first order phase transition at ^ = vr. 



oo. 



§4. Numerical results II 



4.1. Flattening of free energy 

Figure 7 shows F{6) obtained by "direct method" at /? = 3.45. F{0) shows 






1-15 




1-18 




1-20 




1-25 



I ' ' ' ' I ' ' ' ' I ' ' " I ' ' ' ' I " ' ' I ' ' ' ' I ' " ' I 
0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 1 
theta/p i 

Fig. 7. Free energy vs theta at /3 = 3.45 and L = 15, 18, 20, 25. Flattening is clearly seen in 
L = 18,20,25. 

"flattening " at some value Of over L = 18. For 9<0f, F{6) is volume indepen- 
dent. For 9f < 6, F{6) is flat. Of decreases as y — > large. (In L = 25, we can not 
calculate F{0) for > O.Tir, because Z{0) becomes negative.) This flattening phe- 
nomena may correspond to similar behavior found by Schierholz. According to him, 
such flattening is the position of the first order phase transition ( because dF/dO 
is "discontinuous" at the sharp flattening of free energy F{0). Our interpretation 
about the flattening found in our numerical calculations is different. The shift of Of { 
flattening position ) with the size of the volume in our calculations is closely related 
to the statistical error existing in the estimate of Z(0) through fourier series. The 
magnitude of statistical error 5P in the simulation of measurements is essentially 
given by 



6P 



N 



(4-1) 



Statistical error AP{Q) at Q in topological charge distribution P{Q) decreases 
rapidly with Q large in the set method, 



AP{Q) r^6Px P{Q) 



(4-2) 
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where 5P = constant, and then AP{Q) is a rapidly decreasing function of Q. The 
largest contribution of AP{Q) comes from that at Q = 0, 



AP{0) r^5Px P(0). 



(4-3) 



Then Z{9) suffers from the error of P{0) with the magnitude of AP(0) ~ 6P. The 
partition function Z{9) is also a rapidly decreasing function of 9 and then Z{0) 
becomes the order of AP{0) ~ 6P at ~ 0/. In such case, rapidly decreasing 
function Z{0) at > Of can not be estimated safely. 

We will present the results of our simulation. In table III we show the flattening 
position 9f for various lattice sizes L and coupling constants /3, /? = 3.45 and 
P = 3.5. 

Table III. The flattening position 9f for various lattice sizes L. We also show the values of 9f x L. 





L 


Of 




3.45 


15 


0.837r 


I2.4n 




18 


0.77r 


12.67r 




20 


0.67r 


12.l7r 




25 


0.5l7r 


12.87r 


3.5 


20 


0.577r 


11.47r 




30 


O.SStt 


11.47r 



From these results, we observe that the flattening position scales as Of (x L^^ = 
V^^l"^ for each /3, where V = I? . What is the reason for this -law? At these 

/3's, -P(Q) is approximately given by gaussian form 



P(Q) oc exp(-Ky(/5) Q^), 



(4-4) 



and the fourier transform Z{9) of this gaussian form is also gaussian form Z{9) = 
exp(— ay^^). Then free energy F(0^ = — '^lnZ(0) is given by quadratic function 
O'^ . For L = 15 to L = 25, calculated F(G) is almost ^-independent. Since Z(0^ = 
YliQ P{Q)e^^^ suffers from the statistical error 5P, flattening is expected at the 
where the free energy becomes the order of magnitude, 



VF{Of) ~ -ln((5P) 



(4-5) 



and we obtain the relation 

Of ~ y-i/2_ (4.g) 

Of course, 5P may lead to negative value for Z{0) at 9 > Of. In /3 = 3.45, however, 
L = 15, 18, and 20 simulations with = 5 x 10^ iterations all give positive Z{0) 
for all 0(0 < 9 < it). Only L = 25 case leads to negative Z{0) at > O.Tvr. 
We would expect both positive excess fluctuation and negative excess one. So why 
positive excess results were more often observed in /? = 3.45 simulations remains 
as a question. Numerical results for Z{9) and F{0) is shown in the table IV for 
P = 3.45 and L = 20 case. The value of Z{0) in flat region is the same order as 
6P ~ 1/Vn ~ 1.4 X 10-^{N = 5 X 10^). 

Once we found the flattening phenomena in rather large /3 regions, we investi- 
gated the strong coupling region again to see whether it is the phenomenon peculiar 
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to weaker coupling regions. Surprisingly, similar phenomenon is found also in strong 
coupling regions. We show the case of /? = 1.0, which is sufficiently strong; chi-square 
fitting with gaussian form is quite good(Table.I). Fig.l shown in section 1 is the re- 
sult obtained for strong coupling region, /? = 1.0 and L = 10 for 500k iterations. 
Again the position of flattening is consistent with the value estimated from 6P. As is 
seen above, the calculation of the partition function with the (2-8) using the observed 
value of P{Q) directly will be called "direct method". In this method Z{9) receives 
statistical error of the order 5P, and flattening will occur at 9f. On the other hand 
in the set method, the topological charge distribution itself is obtained over quite 
wide range of Q. We first fit it by appropriate analytical function and using the 
fitted function Ps{Q){ "smooth"), we can evaluate Z(9). This latter method will be 
called "fitting method". In strong coupling regions, the indirect method will lead to 
gaussian function for P{Q) and the fourier series will again gives gaussian function 
of 9 for Z{9). It gives smooth 9^ type free energy distribution up to ^ = tt. Namely 
no flattening phenomena will appear. 

Actually in ^(1) LGT, it is known P{Q) is gaussian for all P, and 9^ type free 
energy distribution up to ^ = tt. But if we use the observed P{Q) directly to obtain 
partition function, it will cause flattening phenomena. 

4.2. Strong and weak coupling region 

As is seen in Table.I, In P{Q) can be fitted by gaussian form well up to /? = 
1.5, but beyond that value the fitting becomes worse, and higher power terms are 
necessary. Up to /? = 3.5 quartic function fits rather well, still the most dominant 
power is 2 in this region( see Table. II). In /3 > 4.0 the power series up to the power 
4 is not enough as is seen from Table.II. The value of chi-square become quite bad 
in /? > 4.0. In weaker coupling region (/? > 4.0), we tried another type of function, 
which is series starting from some power 7. 

lnP{Q) = al + Y.al\Q\'r+-' (4-7) 

n=l 

This drastically improves the chi-square fit in 4.0 > /3 > 5.0 (see Table.V. and 
compare with Table.II.). The chi-square value for /? = 4 was 3497 in Table.II, it 
becomes 38.98 in Table.V. Prom Table.V, we understand that leading power 7 is 2 

for (3 < 3.5 and the leading power becomes 1 in /3 > 5, thus we will call this region 
weak coupling region. /3 = 4.0 is "intermediate", since leading power 7 is about 
1.59, which is just medium value of 7 =2(gaussian in strong coupling regions) and 
7=1 (weak coupling regions). 

In weak coupling regions excitation of topological charge is suppressed stronger 
as P becomes larger. For example, 

P(2) _ 1 1 1 
P(l) ~ 132' 3^' 126 

for /3 = 2.0, 4.0, 5.0 respectively. 

On the contrary large topological charge configuration is relatively important 
and gaussian distribution is obtained in strong coupling region. It leads to the 



(4-8) 



14 



M. Imachi, S. Kanou and H. Yoneyama 



partition function, 'i?3(z/, r), third elliptic theta function. Here v = 6/{2tt),t = 
iKv{P)/7T. Note that 6 is extended to complex values. This function clearly shows 
the distribution of partition function zeros leading to the first order phase transition 
at 6 = TT. 

In weak coupling regions large topological charge contribution is suppressed. In 
this region — InP(Q) is given approximately by a7|<3|( with 7 ~ 1.0). This form of 
P{Q) leads to the partition function given by the generating function of Chebyshev 
Polynomials. 

00 

Z{9)= ^ cl^le^^« 
Q=—oo 

_ l-c^ 

~ 1 - 2ccos6' + c2 

In weak coupling regions, the value of c is quite small and cP in this equation can 
be discarded. It leads to the expectation value of topological charge almost given by 
sin^ form. 

i<Q>=l,^/Z^^sine (4-11) 

when |c| « 1. At /3 = 6.0, L = 20, the value of c is c = 1.34 x 10"*^, 2c/V ~ 6.68 x 
10~^. Actually calculated i < Q > shows sine function type behavior over whole 
range of 9 in weak coupling regions. No discrete phase transition is expected at ^ = vr 
in weak coupling regions. Whether there is a continuous phase transition( second 
order one) or not needs further precise investigation such as study of correlation 
length. 

§5. Conclusoins and discussion 

In this paper we have investigated two dimensional CP^ model with ^?-tcrm 
numerically. Topological charge distribution is measured by Monte Carlo method 
based on set method and trial function method. Then it is used to obtain partition 
function Z{9) as a function of theta parameter. 

Topological charge distribution P{Q) is gaussian in strong coupling regions ( 
(3 < 1.5). In /? = 2 ~ 3.5 region, the exponent of P{Q) can not be simply given 
by but additional higher powers are necessary (Table I and II). Still the leading 
contribution comes from term(Table V). 

In weak coupling regions, i.e., (5 > 4.5, the leading power of exponent of P{Q) 
is \Q\. Sharp decrease of P{Q) at larger \Q\ is observed. 

At /3 = 4, which is intermediate of strong and weak coupling region is charac- 
terized by leading power of exponent of P{Q) is given by |(5|^'^^, whose exponent is 
the intermediate value of gaussian P{Q) ) in strong coupling region and \Q\^ in 
weak coupling region. 

Free energy as a function of 9 is investigated. We observed "flattening " of 9- 
distribution of free energy at 9f. In our analysis, the main reason of this flattening 



(4-9) 
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is attributed to the statistical fluctuation 5P of P{Q) at Q = and 5P ^ 0{l/\fN), 
where N is the number of measurements. Volume dependence at each (5 is explained 
by this interpretation. 

Topological charge distribution P{Q) is measured over wide range of Q. But 
direct method to obtain Z{9) through fourier series with direct use of measured 
topological charge distribution P{Q) leads to flattening caused by SP. We tried to 
obtain approximate form of P{Q) by fitting method. The smooth function P{Q)s 
obtained by fitting and used to evaluate Z{9)s has following properties. In strong 
coupling region, the gaussian topological charge distribution leads to Z{9)s with the 
third elliptic theta function ■!?3 and it leads to infinite number of zeros of partition 
function leading to the first order phase transition at = vr in infinite volume limit. 
In weak coupling region, Z{9)s is approximated by the functional form 



which is the generating function of Chebyshev polynomials. Due to this form, the 
expectation value of topological charge 

i<Q>(xsm.9 (5-2) 

No first order phase transition at ^ = tt is expected. 
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Appendix A 

Integration Measure 



In this appendix, we will discuss briefly about the measure of integration. The 
measure is defined as 

dzidz2dzi6{\zi\^ + \z2\'^ + |z3p - 1). (A-1) 

The complex numbers Zi{i = 1,2,3) are given by two real numbers as Zi = 

Xi + iyi- The measure becomes 

dxidyidx'2dy2dzzdz^5{x\ + yl + xl + yl + xl + yl - 1). (A-2) 

When the parameters Xi,yi{i = 1,2,3) are changed to polar coordinates {ki,ai{i = 
1,2,3)) 

Xi = ki cos Oj, (A-3) 
yi = ki sin Qi, (A-4) 

where A;^ > and < Oj < 27r. The measure becomes 

ki k2 ^3 dki dk2 dk^ S{kf + A;| + fc| — 1) dai da2 da^. (A--5) 
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We translate the parameters from ki, ^2, ^3 to 3 dimensional polar coordinates {I, ii, ^2) 

fci = / sinti cos t2, 

^2 = /sintisini2, (A' 7) 

fe3 = /cosi2, (A-8) 

6{P — 1) sin'^ ti cos tidti sint2 cos ^2 dl dai da2 da^, (A-9) 

where / > and < ^1,^2 < f • So the measure which we use in Monte Carlo 
simulation becomes finally 

dm dn dai da2 da^, (A- 10) 

after integration about I is done. The variables m, n are given by 

m = ■^(cos2t2 - 7COs4ti), (A-11) 
8 4 

n= -cos 2*2. (A- 12) 

4 
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Table IV. The free energy and partition function for /3=3.45 and L=20. 
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Table V. The result of chi-square-fitting to \nP{Q) in term of the series starting from Aq + 
E„=i«;ilQr+""' for various 0. 
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